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Abstract 

Previous work in the literature has studied gravitational radiation in black-hole collisions at the 
speed of light. In particular, it had been proved that the perturbative field equations may all 
be reduced to equations in only two independent variables, by virtue of a conformal symmetry 
at each order in perturbation theory. The Green function for the perturbative field equations 
is here analyzed by studying the corresponding second-order hyperbolic operator with variable 
coefficients, instead of using the reduction method from the retarded flat-space Green function 
in four dimensions. After reduction to canonical form of this hyperbolic operator, the integral 
representation of the solution in terms of the Riemann function is obtained. The Riemann function 
solves a characteristic initial-value problem for which analytic formulae leading to the numerical 
solution are derived. 
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I. INTRODUCTION 



In the theory of gravitational radiation, it remains true what R. Sachs stressed in his 
1963 Les Houches lectures, i.e. that we understand the following three main features Q]: 

(i) We can give a description of radiation at large distances from its sources in an asymp- 
totically flat universe; this description is geometrically elegant and sufficiently detailed to 
analyze all conceptual experiments concerning the behaviour of test particles or test ab- 
sorbers in the far field. 

(ii) How the exact theory relates the far field to the near field and the sources or non- 
gravitational fields. 

(iii) We have approximation methods that make it possible to obtain numerical results for 
the amount of radiation emitted in a particular situation, or for scattering cross-sections, 
etc.; these approximation methods do not have a profound geometrical nature, but are very 
important in comparing theory with the experiment, in case the latter succeeds in finding 
observational evidence in favour of gravitational waves. 

Within this framework, it is often desirable to use Green-function methods, since the 
construction of suitable inverses of differential operators lies still at the very heart of many 
profound properties in classical and quantum field theory. For example, the theory of small 
disturbances in local field theory can only be built if suitable invertible operators are consid- 
e red Q. In a in tegral tannic, these „ d to th e g a ug e- fl e ld aed g no st 

operators, respectively ^, 0] . Moreover, the Peierls bracket on the space of physical observ- 
ables, which is a Poisson bracket preserving the invariance under the full infinite-dimensional 
symmetry group of the theory, is obtained from the advanced and retarded Green functions 
of the theory via the supercommutator function 0, [3, U, [J and leads possibly to a deeper 
approach to quantization. Last, but not least, a perturbation approach to classical general 



/ity relies heavily on a careful construction of Green functions of operators of hyper- 
ODDQ and elliptic jlfj, LLH type. In particular, following , we shall be 



concerned with the axisymmetric collision of two black holes travelling at the speed of light, 
each described in the centre-of-mass frame before the collision by an impulsive plane-fronted 
shock wave with energy fi. One then passes to a new frame to which a large Lorentz boost is 
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applied. There the energy v = /ze Q of the incoming shock 1 obeys v » A, where A = /ie~ a 
is the energy of the incoming shock 2 and e a = (P being the usual relativistic param- 

eter). In the^boosted frame, to the future of the strong shock 1, the metric can be expanded 
in the form 



g a b ~ v 2 



8=1 



where 7]^ is the standard notation for the Minkowski metric. The task of solving the 
Einstein field equations becomes then a problem in singular perturbation theory, having to 
find , h^j , ... by solving the linearized field equations at first, second, ... order respectively 
in once that characteristic initial data are given just to the future of the strong shock 1. 
The perturbation series (1.1) is physically relevant because, on boosting back to the centre- 
of-mass frame, it is found to give an accurate description of space-time geometry where 
gravitational radiation propagates at small angles away from the forward symmetry axis 
9 = 0. The news function Co (see appendix), which describes gravitational radiation arriving 
at future null infinity in the centre-of-mass frame, is expected to have the convergent series 
expansion 

oo 

co(f,0) = ^a 2n (f//i)(sin0) 2n , (1.2) 

n=0 

with f a suitable retarded time coordinate, and \i the energy of each incoming black hole in 
the centre-of-mass frame. In a very useful analytic expression of a 2 (f/fi) was derived, 
exploiting the property that perturbative field equations may all be reduced to equations 
in only two independent variables, by virtue of a remarkable conformal symmetry at each 
order in perturbation theory. The Green function for perturbative field equations was then 
found by reduction from the retarded flat-space Green function in four dimensions. 

However, a direct approach to the evaluation of Green functions appears both desirable 
and helpful in general, and it has been our aim to pursue such a line of investigation. For 
this purpose, following hereafter our work in reduction to two dimensions with the 
associated hyperbolic operator is studied again in section 2. Section 3 performs reduction to 
canonical form with the associated Riemann function. Equations for the Goursat problem 
obeyed by the Riemann function are derived in section 4, while the corresponding numerical 
algorithm is discussed in section 5. Some backgound material is described in the appendix. 
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II. REDUCTION TO TWO DIMENSIONS AND THE ASSOCIATED OPERATOR 

As is well known from the work in Q and Q, the field equations for the first-order 
correction in the expansion (1.1) are particular cases of the general system given by the 
flat-space wave equation (here u = -A=(z + t),v = -4^{z — £)) 

_ 2 d 2 ip + 1 d ( dif>\ + 1 d 2 ip _ Q 
owat> p dp \ dp J p 2 d<p 2 

supplemented by the boundary condition 

^(u = 0) = e m V n /[8 logM - y/2v), (2.2) 

/(ar) = Va; < 0. (2.3) 
Moreover, if) should be of the form e tm ^ p~ n \(q 1 r) for w > 0, where 

g = M p" 2 , (2.4) 



r = 8 log(i/p) - V2v. (2.5) 

For the homogeneous wave equation (2.1) there is no advantage in eliminating p and (f> from 
the differential equation. However, the higher-order metric perturbations turn out to obey 
inhomogeneous flat-space wave equations of the form 

□ V = S, (2.6) 

where S is a source term equal to e irn( ^p~^ n+2 ^H(q,r). This leads to the following equation 

for x = e~ im ^p n ifj: 

An,nX(g,r) = H(q,r), (2.7) 
where C™ n is an hyperbolic operator in the independent variables g and r, and takes the 



form 



^2 o2 Q2 

C m , n = -(2 V / 2 + 32g) — - + 4g 2 — + 64— 

agar ag z ar z 

+ 4(n + l)g— - I6n— + n 2 - m 2 . (2.8) 
ag or 
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The proof of hyperbolicity of Criyn, with the associated normal hyperbolic form, can be 
found in section 3 of 3], and in p. The advantage of studying Eq. (2.7) is twofold: to 
evaluate the solution at some space-time point one has simply to integrate the product of 
H and the Green function G m , n of C m ^ n : 

X(V, r ) = J G min (q, r; q , r )H(q , r )dq dr , (2.9) 



and the resulting numerical calculation of the solution is now feasible 
If one defines the variables 



X = log(g) + ^, y = log(g)-^ (2.10) 

the operator C m , n is turned into 

d 2 d 1 / d 2 d 2 \ 
T = 16—— + 8n— + n 2 -m 2 ^ e ~ (x+y)/2 I — — I (2 11) 

The operator T m n is the 'sum' of an elliptic operator in the Y variable and a two-dimensional 
wave operator 'weighted' with the exponential e~( x+Y ^ 2 , which is the main source of tech- 
nical complications in these variables. 



III. REDUCTION TO CANONICAL FORM AND THE RIEMANN FUNCTION 

It is therefore more convenient, in our general analysis, to reduce first Eq. (2.7) to 
canonical form, and then find an integral representation of the solution. Reduction to 
canonical form means that new coordinates x = x(q,r) andv = y(q,r) are introduced such 
that the coefficients of and vanish. As is shown in 2,0], this is achieved if 



dx dy 

dr dr ' ^ 



dx _ l + 8gy/2+ ^1 + 16^^2 

dq~ 2V2q 2 ' { } 



dy l + 8qy/2-y/l + 16qV2 



dq 2y/2q 2 
The resulting formulae are considerably simplified if one defines 



t = Jl + 16qV2 = t(x,y). (3.4) 



The dependence of t on x and y is obtained implicitly by solving the system 

x = r + log ^— — j - - 4, (3.5) 



This leads to the equation 



*= r+to «^) + (Fnr 4 (3 - 6) 



2t (x-y) 

log (^Ti)-(^i) = ^' (3 - 7) 



which can be cast in the form 

(^ e <^>=e^. ( 3. 8 ) 



+ 

This suggests defining 



(t+iy 

so that one first has to solve the transcendental equation 



w = -. — -, (3.9) 



\ w ~ i ) \ x ~y) 

we as = , (3.10) 
to obtain w = w(x — y), from which one gets 

(1-w) 

On denoting by g(w) the left-hand side of Eq. (3.10), one finds that, in the plane (w, g(w)), 
the right-hand side of Eq. (3.10) is a line parallel to the w-axis, which intersects g(w) at no 
more than one point for each value of x — y. For example, when w — 1, g(w) intersects the 
line taking the constant value 1, for which x — y = 0. The function 

g : w —> g{w) = we 2 ™ 

is asymmetric and has the limiting behaviour described by 

lim g{w) = — oo, lim g(w) = 0, (3-12) 

w— >0 _ w— >0+ 



lim = 0, lim g{w) = oo. (3.13) 
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Thus, in the lower half-plane, g has an horizontal asymptote given by the w-axis, and a 
vertical asymptote given by the line w — 0, while it has no asymptotes in the upper half- 
plane, since 

hm = oo 

w— »oo W 

in addition to (3.13). The first derivative of g reads as 



[w 



l) 2 (w *-p 



2w 



(3.14) 



One therefore has g'(w) > for all w > 0, and g'(w) < for all w G (— oo, 0) — {—1}, and 
g is monotonically decreasing for negative w and monotonically increasing for positive w. 
The point w — —1, at which g'(w) vanishes, is neither a maximum nor a minimum point, 
because 

(3.15) 



l l i w . 

1 h 1 H e 2™ 



1 3 3 

+ 



3 3 w\ (-w 2 -i) 
e 2 ™ . 



8w 5 4u> 4 8w 3 8w 4 
These formulae imply that g"(—l) = but g"'(—l) = —1^0, and hence u> 
flex of g(w) (see Fig. 1). 

In the (x, y) variables, the operator C m ^ n therefore reads 

d 2 d d 



where, exploiting the formulae 



one finds 



dx 6AV2 



dq (t-1) 2 ' 
dy 64V2 



dq (t + 1) 2 ' 



f(x,y) 



= 256 



2t 2 (t 2 + 1) 



dx dy 
dq dq 



+ 128 



{t-mt+if 



(3.16) 
-1 yields a 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



dx 

g(x, y) = 4(n + l)q — 16n = 16 



1 + 



2(w+ 1) 
(t-1) 



(3.21) 
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FIG. 1: Plot of the function g : w — > <?(ii>) = we 



(u> 2 -l)/2u; 



/i(x, y) = 4m + 16n = 16 



2(w+ 1) 
+ 



The resulting canonical form of Eq. (2.7) is 

1- n\ t. ii t 

dx 



where 



( d 2 d d 

I + a(x, y)— + 6(x, y)— + c(x, y) ) y) 



\dxdy 
H(x,y) 



dy 



a{x,y) 

b(x,y) 
c(x,y) = 
H(x,y) = 



g(x,y) _ 


1 (1- 


-t)(t + l) 2 (2n + l + t) 


f(x,y) 


16 


(t 4 + At 2 - 1) 


h(x,y) _ 


1 (t + 


l)(t - l) 2 (2n+ 1 - t) 


ffay) 


16 


(t 4 + 4t 2 - 1) 


i i 
n — m 


(m 2 


-n 2 ) (t-l) 2 (t+l) 2 



H(x,y) 
f(x,y) 



256 (t 4 + 4t 2 -l) 

H(x,y) (t-l) 2 (t + l) 2 
256 (t 4 + 4t 2 -l) ' 



(3.22) 



(3.23) 

(3.24) 
(3.25) 
(3.26) 
(3.27) 
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Note that a(-t) = b(t),b(-t) = a{t),c(-t) = c(t),H(-t) = H(t). 

For an hyperbolic equation in the form (3.23), we can use the Riemann integral represen- 
tation of the solution. For this purpose, recall from Q| that, on denoting by the adjoint 
of the operator L in (3.23), which acts according to 

L ] [x] = Xxy ~ (ax)x ~ My + ex, (3.28) 

one has to find a 'function' R(x, y; £, rj) (actually a kernel) subject to the following conditions 
((£,77) being the coordinates of a point P such that characteristics through it intersect a 
curve C at points A and B, AP being a segment with constant y, and BP being a segment 
with constant Xj cLS IS shown in Fig. 2): 

(i) As a function of x and y, R satisfies the adjoint equation 

L\ X JR] = 0, (3.29) 

(ii) R x = bR on AP, i.e. 

R x (x, y; C> v) = K x , v)R(%> y\ C> »?) on 2/ = v, (3.30) 
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and R y = aR on BP, i.e. 



Ry(x, v\ 6, v) = a(f , y; £, r?) on x = £, 

(iii) i? equals 1 at P, i.e. 

It is then possible to express the solution of Eq. (3.23) in the form 



X(P) = 1 -[x(A)R(A) + X (B)R(B)]+ [ 
1 Jab 



^Xx + (bR - ^R x ) \ 



dx 



T^Xy + (aR - ^Ry ) \ 




(3.31) 



(3.32) 



dy)+ I I R(x,y;£,r])H(x,y)dxdy, (3.33) 



where Q is a domain with boundary. 

Note that Eqs. (3.30) and (3.31) are ordinary differential equations for the Riemann 
function R(x, y; £, if) along the characteristics parallel to the coordinate axes. By virtue of 
(3.32), their integration yields 

/X 
b(\,r))d\, (3.34) 



R(£,y,£,v) = a(\,£)d\, (3.35) 

which are the values of R along the characteristics through P. Equation (3.33) yields 
instead the solution of Eq. (3.23) for arbitrary initial values given along an arbitrary non- 
characteristic curve C, by means of a solution R of the adjoint equation (3.29) which depends 
on x, y and two parameters £, rj. Unlike x-, the Riemann function R solves a characteristic 
initial-value problem. 



IV. GOURSAT PROBLEM FOR THE RIEMANN FUNCTION 

By fully exploiting the reduction to canonical form of Eq. (2.7) we have considered novel 
features with respect to the analysis in |?[ E|, because the Riemann formula (3.33) also 
contains the integral along the piece of curve C from A to B, and the term ^[x(A)R(A) + 
x{B)R(B)]. This representation of the solution might be more appropriate for the numerical 
purposes considered in [8], but the task of finding the Riemann function R remains extremely 
difficult. One can however use approximate methods for solving Eq. (3.29). For this purpose, 
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we first point out that, by virtue of Eq. (3.28), Eq. (3.29) is a canonical hyperbolic equation 
of the form 



where 

A = -a, (4.2) 

B = -6, (4.3) 

C = c-a x -b y . (4.4) 

Thus, on defining 

U = R, (4.5) 

V = R X + BR, (4.6) 
the equation (4.1) for the Riemann function is equivalent to the hyperbolic canonical system 

Q 

U x = h{x,y)U + h{x,y)V, (4.7) 

V y = g 1 (x,y)U + g 2 (x,y)V, (4.8) 

where 

fi = ~B = b, (4.9) 

h = 1, (4.10) 

g 1 = AB - C + B y = ab - c + a xi (4.11) 

g 2 = -A = a. (4.12) 
For the system described by Eqs. (4.7) and (4.8) with boundary data (3.34) and (3.35) 



an existence and uniqueness theorem holds (see 



13j for the Lipschitz conditions on bound- 



ary data), and we can therefore exploit the finite differences method to find approximate 
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solutions for the Riemann function R(x,y,£,r)), and eventually x{P) with the help of the 
integral representation (3.33). 



V. CONCLUDING REMARKS 

I — I 

The inverses of hyperbolic operators |l4j ] and the Cauchy problem for hyperbolic equations 
with polynomial coefficients 15[ have always been the object of intensive investigation in 
the mathematical literature. We have here considered the application of such issues to 
axisymmetric black hole collisions at the speed of light, relying on the work in 
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9[. We have pointed out that, for the inhomogeneous equations (2.7) occurring in the 
perturbative analysis, the task of inverting the operator (2.8) can be accomplished with the 
help of the Riemann integral representation (3.33), after solving Eq. (4.1) for the Riemann 
function. One has then to solve a characteristic initial-value problem for a homogeneous 
hyperbolic equation in canonical form in two independent variables, for which we have 
developed formulae to be used for the numerical solution with the help of a finite differences 
scheme. For this purpose one studies the canonical system (cf (4.7) and (4.8)) 

U x = F(x, y, U, V), (5.1) 



V y = G(x,y,U,V), (5.2) 

in the rectangle 1Z = {x, y : x G [x , xq + a], y G [yo, yo + b]} with known values of U on the 
vertical side AD where x = xq, and known values of V on the horizontal side AB where 
y = yo. The segments AB and AD are then divided into m and n equal parts, respectively. 
On setting — = h and - = k, the original differential equations become equations relating 
values of U and V at three intersection points of the resulting lattice, i.e. 

U(P M ) - U{P rs 



h 



F, (5.3) 



k 

It is now convenient to set U rs = U(P rs ), V rs = V(P rs ), so that these equations read as 

U r ,s+i — U rs + hF{P rs , U rs , V rs ), (5-5) 
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V r+X ,s = V rs + kG(P rs , U rs , V r8 ). (5.6) 

Thus, if both U and V are known at P rs , one can evaluate U at P r , s +i and V at P r +i, s - 
The evaluation at subsequent intersection points of the lattice goes on along horizontal or 
vertical segments. In the former case, the resulting algorithm is 

s-1 

U rs = U r0 + hY^ F(P ri , U ri , V ri ), (5.7) 
i=i 

V rs — Vr~i,s + kG(P r _i tS , J7 r _x, s , V r - l, s ), (5-8) 
while in the latter case one obtains the algorithm expressed by the equations 

r- 1 

V rs = V 0s + J2 G ( p *s, U is , Vu), (5.9) 
1=1 

U rs = U r ^ x + hF(P r>s _ x , U rta - U (5-10) 

Stability of such solutions is closely linked with the geometry of the associated characteristics, 
and the criteria to be fulfilled are studied in section 13.2 of [l(| (stability depends crucially 
on whether or not | < 1). 

To sum up, one solves numerically Eq. (3.10) for w = w(x, y) = w(x — y), from which one 
gets t(x — y) with the help of (3.11), which is a fractional linear transformation. This yields 
a, b, c and H as functions of (x, y) according to (3.24)-(3.27), and hence A, B and C in the 
equation for the Riemann function are obtained according to (4.2)-(4.4), where derivatives 
with respect to x and y are evaluated numerically. Eventually, the system given by (4.7) 
and (4.8) is solved according to the finite-differences scheme of the present section, with 

F = fiU + f 2 V = f x R + f 2 (R x + BR), (5.11) 



G = gi U + g 2 V = gi R + g 2 (R x + BR). (5.12) 

Once the Riemann function R = U is obtained with the desired accuracy, numerical evalu- 
ation of the integral (3.33) yields x(P)i an d x(<?5 r ) is obtained upon using Eqs. (3.5) and 
(3.6) for the characteristic coordinates. 
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Our steps are conceptually desirable since they rely on well established techniques for 
the solution of hyperbolic equations in two independent variables Q], and provide a 
viable alternative to the numerical analysis performed in because all functions should 
be evaluated numerically. Our method is not obviously more powerful than the one used 
in , but is well suited for a systematic and lengthy numerical analysis, while 

its analytic side provides an interesting alternative for the evaluation of Green functions 
both in black hole physics and in other problems where hyperbolic operators with variable 
coefficients might occur. This task remains very important because a strong production 



of gravitational radiation is mainly expected in the extreme events studied in 
and which motivated our paper. Any viable way of looking at mathematical and numerical 
aspects of the problem is therefore of physical interest for research planned in the years to 
come 111 



APPENDIX A: THE CHARACTERISTIC INITIAL- VALUE PROBLEM 

In our expository article, we find it appropriate to include some background material, 
following, for example, the presentation in section IV of We therefore consider some four- 
dimensional region of space-time and choose in it a set of null hypersurfaces u = constant; 
the corresponding ray congruence with tangent vector k a = u^ a is assumed to have expansion 
p = \k a . a 7^ 0, which can always be arranged in a space-time patch, whereas outside of some 
patch the rays start to cross and hence our construction breaks down globally. On completing 
the k a direction to a quasi-normal tetrad (k, m, t), one finds the following split of the vacuum 
Einstein equations with Einstein tensor Gab- 
Main equations (6 equations) 

k a G ab = 0, G ab t a t b = 0, (Al) 

trivial equation 

G ab t a t = 0, (A2) 

and 3 supplementary conditions 

G ab m a t b = 0, G ab m a m b = 0, (A3) 

where a single complex equation has been counted as two real equations. Remarkably, if 
the main equations hold everywhere, then the trivial equation holds everywhere and the 
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supplementary conditions hold everywhere if they hold at one point on each ray. The 
fulfillment of the trivial equation is proved by writing, from the vacuum Einstein equations, 
that 

G\ a = 0, (A4) 
and then exploiting the main equations (Al) jointly with the split of k a . b as given in Q]: 

K;b = zt a tb + crt a t b + Qt a k b + (k a t b + c.c. 

+ £k a k b . (A5) 



Hence one gets 



= k a G\. h = -k a , b G ab = -2pG ab t a t. (A6) 



By hypothesis the expansion p does not vanish, so that the trivial equation is, indeed, 
identically satisfied. The fulfillment of (A3) everywhere is proved along similar lines. Thus, 
one can again integrate the main equations (Al) first and worry about the supplementary 
conditions (A3) later. 

Choose now the coordinate x° as the retarded time: x° = u. Let r = x 1 be a luminosity 
distance along the rays; let x a (with a = 2, 3) be any other pair of coordinates constant 
along the rays. The line element in these coordinates takes therefore the form (no confusion 
should arise with the j3 of section 1) 

ds 2 = W du 2 + 2e 2/3 dudr - r 2 h ai (dx a - U a duj (dx 1 - W duj , ( A7) 

where W, /3, h ay , U a depend on the x a coordinates. Since r is a luminosity distance, the 
determinant of h Vfl is independent of r. Bearing in mind that the luminosity distance 
is defined only up to a factor constant along each ray, one can demand without loss of 
generality that (here 9 = x 2 , <fi = x 3 ) 

2h flu dx> 1 dx u = (e 27 + e 25 ) d9 2 + 4 sinh( 7 - 5)d9d(j> sin 9 

+ sin 2 fl(V 27 + e- 25 )^ 2 . (A8) 

The metric corresponding to the line element (A7) contains only six unknown functions of 
four variables, and our coordinate system is 'rigid' enough for our purposes [l|. 

One can either analyze the field in the neighbourhood of some point, or the field near 
infinity in an asymptotically flat space-time. Indeed, if in Minkowski space-time one uses a 
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retarded time u = t — r and spherical coordinates r, 6, 4> one finds for the line element 

d s 2 = du 2 + Idudr - r 2 (d6 2 + sin 2 9 dcp 2 ^ , (A9) 
hence one is led to require that, if asymptotic flatness holds, 

lim W =1, lim (rU a ) = lim f3 = lim 5= lim 7 = 0, (A10) 



r— >oo r— *oo r—*oo r— >oo r— >oo 



where all limits are taken as r approaches infinity with u, 9, <p fixed. The second requirement 
in (A10), i.e. that all quantities of interest admit a power-series expansion in -, e.g. 

(1 _,) ( i + , T)/2 = *M + *M + ..., ( Aii) 

is indeed restrictive. Such a requirement can be drastically weakened but not fully elimi- 
nated; moreover, it is closely related to an outgoing radiation condition of the Sommerfeld 
type. It should be stressed that all these requirements no longer hold when r becomes small 
to the extent that rays start to cross each other. 

In the axially- and reflection-symmetric case considered by Bondi et al. jisjl , one has 

5 = 7, ^ 3 = 0, ^ = 0, (A12) 

and the ^-direction is a hyper surf ace-orthogonal Killing direction. The line element acquires 
the simpler form 



Ve 213 

ds 2 = du 2 + 2e 2/3 dudr - r 2 

r 



e 2 \d6 -U du) 2 + e" 27 sin 2 9 d<f) 2 , (A13) 



where the peculiar form of the first coefficient is chosen to simplify the resulting calculations. 

Interestingly, two main equations are found to be identically satisfied by virtue of axial 
symmetry, whereas the other four turn out to be linear combinations of [18( 

Ru = ~(Pi-\ni) =0, (Ai4) 



-2r 2 R 



12 



2r 2 - 712 + 2 7l 7 2 - 27! cot 9 + (3 12 - 2 



0, (A15) 



r 2 U 12 



ArU 2 - r 2 Ui cot 9 - ArU cot 9 



+ 2e 2 ^ 



P22 + 131 - 1 - (3 72 - fa) cot 9 - 722 + 272(72 - #2) = 0, (A16) 



16 



-r 2 R\e w = 2r(r 7 )oi + (1 - r7i)Vi - {ry u + 7i)V - r(l - ry 1 )U 2 
- r 2 (cot 9 - 7 2 )C/i + r(2r7i2 + 2 7a + ryi cot - 3 cot 6")Z7 
+ e 2 ^-^ [ - 1 - (3 72 - 2/? 2 ) cot - 722 + 2 72 (7 2 - /? 2 )j = 0. (A17) 

Equations (A14)-(A16) are called hypersurface equations because they contain no u deriva- 
tives, while Eq. (A17) is called the standard equation. 

Now if 7 is given for one value of u, Eq. (A14) and the boundary conditions (A10) 
determine (3 uniquely. Next Eq. (A15) and the boundary conditions determine U up to 
a function of integration — QN(u, 9) that can be added to r 4 e 2(7_/3) C/i. Equation (A16) 
determines V up to the additive function —2M(u,9); last, Eq. (A17) determines 70 up to 
an additive function co< -"'^ . One can then differentiate Eqs. (A14)-(A17) with respect to 
u and repeat the whole procedure. To sum up, given 7 at one moment the main equations 
determine the future or past up to the three integration functions just mentioned. In the 
general case, the results are completely similar. One has to assign at one value of u the two 
functions 7 and 5. The future is then determined up to five integration functions: a term 
—2 MvM>jO to be added to W; two functions N a (u, 9, 0) which occur in the r~ 3 term for U ; 
and two 'news functions' cq(u, 6, 4>), where the complex function c is given by Eq. (All). 

As far as the supplementary conditions are concerned, the lemma just given makes it clear 
that they should only involve the functions M, N and Co, while a long calculation yields 

I, . f. 



M = -\c \ 2 + ^{sm6)- i Re{ V 



l/sin#)V(c sin 2 #) }, (A18) 



3(N 2 + i sin 6N 3 ) = - VM - [4c cot 9 + (Vc) + 3cV]c , (A19) 

where V = J| + i(sin 6 1 ) -1 ^. The desired M and A^" can be determined once that c(u, 9, <fi) 
and some initial values are given. In the axially symmetric case, Eqs. (A18) and (A19) take 



the simpler form [1 



M = -eg + ^(c 22 + 3c 2 cot 9 - 2c) , (A20) 



-3A" = M 2 + 3c c 02 + 4cc cot 9 + c c 2 . (A21) 

The functions 7 and 5 given on the initial hypersurface u = constant, jointly with the 
two news functions Co given at r = 00 describe the two transverse degrees of freedom. 
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Moreover, one should specify M and N at the initial (or final) retarded time, and these three 
functions of two variables must be related to the longitudinal-timelike degrees of freedom 
of the gravitational field. In the characteristic value problem for general relativity, the 
independent data appear therefore in a very explicit form 
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